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We calculate the photo-ionization cross-section from the ground state of the helium atom, us- 
ing the complex rotation method and diagonalization of sparse matrices. This produces directly 
the positions and widths of the doubly excited P° resonances together with the photo-ionization 
cross-section. Our calculations up to the N — 9 threshold are in perfect agreement with recent 
experimental data and show the transition from a regular structure at low energy to a chaotic one 
at high energy, where various resonances strongly overlap. 



The Helium atom is one of the prototype of atomic systems whose classical dynamics is mainly chaotic and during 



the past thirty years, it has been matter of numerous studies from both theoretical [|l|-|8| and experimental |9|,[10 12 
i points of view. But, unlike other atomic systems like for instance the hydrogen atom in magnetic field, the effects 
pL^ ' of chaos are not very well understood and more profound studies are needed. This requires the resolution of the full 
quantum problem. Especially, one has to take into account all the degrees of freedom and all correlations between 
the two electrons, as well as the auto-ionizing character of the doubly excited states. 

In this letter, we present numerical calculations of the cross-section of the one-photon photoionization from the 
ground state of the helium atom and compare them with the recently obtained high-resolution spectra of 1 P° doubly 
excited states. The agreement with the most recent experimental data from references [ p~0|JT2) - up to the N = 9 
ionization threshold, less than 1 eV from the double ionization threshold and corresponding to 64 open channels - is 
i excellent for the whole energy range, proving the high efficiency of the method. Predictions for better experimental 
resolutions are also given. We also show that at low energy the resonances can be classified with respect to the 
Herrick's (N,K,T) approximate quantum numbers ((N,K) n Lin's simplified notation || will be used hereafter). 
At high energy, this classification progressively breaks down. Eventually, above the N = 7 threshold, the various 
£h ■ resonances strongly overlap: the mean energy spacing between consecutive resonances becomes smaller than their 
typical width. Oscillations in the photo-ionization cross-section can then no longer be associated with individual 
resonances: random-like fluctuations - known as Ericson fluctuations - should be observed in the cross-section. 
O ' The quantum Hamiltonian in atomic units (h = m e - = Aireo = e 2 = 1) is given by: 
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where Pj = —iKVi is the momentum operator of electron i, ri its distance to the nucleus and 7*12 the inter-electronic 
distance. All spin-orbits, relativistic and QED effects (at most of the order of a fraction of meV) are neglected, which 
is consistent with the experimental resolution (of the order of 1 meV). The corrections due the finite mass of the 
nucleus are taken into account by using the effective values for the double-ionization threshold 1^ and the Rydberg 
constant Rn e given in the reference JlO) , namely 1^ = 79.003 eV and Rn e — 13.6038 eV. 

Using the rotational invariance of the Hamiltonian, the angular dependency of a wavefunction can be factorized as 
follows Q: 

L 

*LM= VwiM^T^frV,*) (2) 

T=-L 

where (-0, 0, <fi) are Euler angles defining the transformation from the laboratory frame to a molecular-like frame whose 
z axis is the inter-electronic axis. \T\ is then the A (E, n...) quantum number in a molecule. The T>j^ T are the 
wavefunctions of the rigid rotor and reduce to the usual spherical harmonics for T = 0. Finally (x, y, z) are the 
perimetric coordinates, symmetric combinations of 7*1, r 2 and ri2: 

7*12 

(3) 
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For each pair of good quantum numbers (L,M), we obtain an effective Hamiltonian acting on the different $t's 
(coupled by Coriolis-like terms). The two remaining discrete symmetries - parity and exchange between the two 
electrons - are exactly taken into account by adding constraints on the <I>t's Isl. 

As stated before, above the first ionization threshold, all states become resonances because of the coupling with 
the continua (autoionizing states). Using the complex rotation method |5|,|l3|], we obtain these resonances as complex 
eigenvalues of a complex Hamiltonian H(9), which is obtained by the replacements — > r^e* 6 and Pi — * Pie~ 10 , 
where 9 is a real parameter. The fundamental properties of the spectrum of H (9) are the following : the continua 
of H are rotated by an angle 29 in the lower complex half-plane around their branching point. Each other complex 
eigenvalue Ei = — iTi/2 lies in the lower half-plane and coincides with a resonance of H with energy ej and width 
IV These quantities are independent of 9 provided that the complex eigenvalue has been uncovered by the rotated 
continua. The bound states, which are resonances with zero width, stay on the real axis. This method also allows to 
compute quantities of physical interest, like photo-ionization cross-section, probability densities or expectation values 
of operators (e.g. cos #12), enlighting the contribution of a given resonance to them. For instance, the cross-section is 
given by ||| : 

^) = ^EPf¥ (4) 
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T is the dipole operator, fku is the photon energy, \g) is the ground state (of energy E g ). (Ei$\ is the transpose of the 
eigenvector \Eie) of H(6) for the eigenvalue Em (i.e. the complex conjugate of (Ei$\). R(9) is the rotation operator, 
essential to obtain the right (complex) oscillator strength. 

In the preceding formula, each eigenvalue (resonance or continuum) contributes to the cross-section at energy 
fko + Eg, with a Fano profile centered at energy KeEi, of width — 2Im.Ej whose q parameter is given by |f3( ]: 

= Re(E^\R(9)T\g) 
q lm(E ie \R(9)T\g) 

Thus, the Fano-g parameter of one resonance is directly and unambiguously obtained from its associated eigenvector, 
which is much more efficient than any fitting procedure, especially above the N = 6 threshold where the different 
series strongly overlap (see fig. [I]) . 

For an efficient numerical resolution, the effective Hamiltonian is expanded in the product of three Sturmian-like 
basis \n x ) ® \n y ) ® \n z ), one for each perimetric coordinate. The basis states have the following expression: 

with (u\n) = 4> n (u) = y/a~^L n (a u u)e~ a "2 (6) 

where n x , y ,z a^e n on negative integers, a x ,y,z are rea l positive parameters (the scaling parameters) and L n the n th 
Laguerre polynomial. This non-orthogonal basis is associated with a representation of the dynamical group SO(2,l), 
which gives rise to selection rules. The matrix representation of the effective Hamiltonian in this basis is thus sparse 
and banded, and the matrix elements are analytically known. Let us emphasize that this approach is "exact" for the 
non-relativistic He atom and similar to the one used in reference ||. 

For obvious reasons, the basis has to be truncated, the prescription being n x + n y + n z < N max (we used up to 
JV max = 58). The different scaling parameters are related by a — a x — 2a y = 2a z , which increases the sparsity 
of the matrices and gives the correct decrease for r± and r 2 going to infinity. The matrices are diagonalized with 
the Lanczos algorithm, which is a highly efficient iterative method to obtain few eigenvalues of huge matrices in a 



short CPU time 14 1. Convergence of the results are checked with systematic changes of a and 9. We have thus 
computed few hundred 1 P° states, which are the only ones populated in one photon experiment starting from the 
Helium ground state ( 1 S e ). The resulting cross-section from below the N = 2 up to above the N = 8 threshold - 
the highest energy where experimental spectra are available - is shown in figure |, convoluted with a lorentzian at 
the experimental resolution (2 meV for TV = 2, 3 and 4 meV for N — 4,5, 6) or at a slightly better resolution (1 meV 
above the N = 6 threshold). The agreement with the figures from references @,^2[ is excellent, emphasizing the 
efficiency of our calculations. The theoretical positions, linewidths and Fano-g parameters are in good agreement with 
previous works [g 10|. 

Below the N — 2 (resp. N = 3), three (resp. five) different series are clearly distinguishable, either by their widths 
(see fig. |l|) or by the expectation value of cos6> 12 , as shown in Fig. ||[ where the real part of — N cos9i 2 is plotted (the 
imaginary part is at least ten times smaller) versus the effective principal quantum number n e g of the outer electron 
measured from the N th threshold, proving thus the validity of Herrick classification in these energy ranges. Still, the 
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chaotic aspect of the Helium atom is already observable in the fluctuations of the smallest widths (see Fig. P (and 
also in the Fano-g parameters), which will be amplified at higher energies. 

Below the N = 5 and N = 6 thresholds, irregularities due to the interaction with the 6,4g (resp. 7,57) state from 
the upper series are visible, in perfect agreement with the experimental observation. Below the N = 7, N — 8 and 
N = 9 thresholds, the computed cross-section - represented at a better resolution - reproduces very well the various 
overlapping series, with an increasing number of perturbers coming from higher series. Furthermore, we show new 
peaks that are not yet experimentally resolved - such as the members of the 9, 7 n series - but whose observation 
could be possible with a (slight) increase of the experimental resolution and signal to noise ratio. In this energy 
range, the various series are so strongly coupled and overlapping that the approximate Herrick classification breaks 
down [ p^lp^ ] , giving rise to an irregular spectrum, showed by Fig. [I] where no general trend can be easily recognized in 
the widths of the various resonances. This irregularity is the quantum manifestation of the chaotic classical dynamics. 
In this regime, the photo-ionization cross-section results from the superposition of various overlapping Fano profiles, 
eventually leading to random-like fluctuations in the cross-section known as Ericson fluctuations |17| . Predicted 
around the N = 30 threshold in the 1-dimensional helium atom fl8[] , this irregular regime takes place at much lower 
energy in the real helium atom, because of the increased density of states. The ratio between the linewidth T and 
the local mean level spacing s is displayed in Fig. || for the N = 4, N = 6 and N — 8 thresholds. We clearly see 
that for N = 8, a vast majority of resonances lie above the line T/s — 1, corresponding to the overlapping resonances 
regime. The published experimental results Jl2|, seem to show irregular fluctuations, the first steps towards the 
Ericson fluctuations. 

In conclusion, our results are, as far as we know, the ab initio calculations for the double excited 1 P° states 
of the helium atom at the highest energy ever done. They are in excellent agreement with the presently available 
experimental data. Importantly, they show that the strongly irregular regime where various resonances overlap leading 
to Ericson fluctuations in the photo-ionization cross-section is almost reached experimentally, which opens the way 
to their experimental observation and more generally to a new generation of experiments probing the chaotic aspects 
of the helium atom. 
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FIG. 1. Positions (in eV above the ground state) and widths (logarithmic scale) of 1 P° resonances of the Helium atom. 
The upper plot displays the states below the N = 2 and iV = 3 thresholds, where the various series can be distinguished 
without ambiguity from their widths, in agreement with Herrick classification of doubly excited states. In the lower plot, 
displaying the states below the N = 6, 7, 8 thresholds, the various series are strongly coupled and overlapping, which is a 
quantum manifestation of classical chaos in this system. 
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FIG. 2. Calculated photo-ionization cross-section of the helium atom from the N — 2 (upper plot) to the N = 8 and 9 (lower 
plot) series. The raw spectrum has been convoluted by a lorentzian of width 2 meV for JV = 2, 3 and 4 meV for iV = 4, 5, 6 
(equal to the experimental resolution) and 1 meV for N = 7, 8, 9. The calculated cross-section is in excellent agreement with 
the experimental results of references |l^,^2), displaying for N = 7, 8, 9 new peaks, not yet experimentally observed. At the 
highest energies, the various series overlap strongly, leading to irregular fluctuations of the cross-section and breakdown of the 
classification. Only the fluctuating part of the cross-section is here represented, the smooth background being subtracted. 
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FIG. 3. Real part of the expectation values of —N cos 612 for the various resonances below the N = 2 and N = 3 thresholds. 
As expected from the Herrick classification scheme of doubly excited states, the value is almost constant across a series, although 
it slightly differs from the predicted value — (N — 1) < K < (N — 1). At higher energy, this classification breaks down. 
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FIG. 4. Ratio between the widths V of the various 1 P° resonances of the helium atom and the local mean level spacing s, 
displayed for the N — 4, N — 6 and N = 8 thresholds (from left to right). The transition between the regime of well separated 
resonances to the strong overlapping resonances regime is observed. For higher thresholds, the number of resonances lying 
above the T = s line will increase, leading to the observation of Ericson fluctuations in the photo-ionization cross-section. 
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